Spatially explicit effective reproduction numbers from incidence and mobility data.

Current methods for near real-time estimation of effective reproduction numbers from surveillance data overlook mobility fluxes of infectors and susceptible individuals within a spatially connected network (the metapopulation). Exchanges of infections among different communities may thus be misrepresented unless explicitly measured and accounted for in the renewal equations. Here, we first derive the equations that include spatially explicit effective reproduction numbers, ℛk(t), in an arbitrary community k. These equations embed a suitable connection matrix blending mobility among connected communities and mobility-related containment measures. Then, we propose a tool to estimate, in a Bayesian framework involving particle filtering, the values of ℛk(t) maximizing a suitable likelihood function reproducing observed patterns of infections in space and time. We validate our tools against synthetic data and apply them to real COVID-19 epidemiological records in a severely affected and carefully monitored Italian region. Differences arising between connected and disconnected reproduction numbers (the latter being calculated with existing methods, to which our formulation reduces by setting mobility to zero) suggest that current standards may be improved in their estimation of disease transmission over time.

Assessing the impact of non-pharmaceutical interventions on SARS-CoV-2 transmission in Switzerland.

Following the rapid dissemination of COVID-19 cases in Switzerland, large-scale non-pharmaceutical interventions (NPIs) were implemented by the cantons and the federal government between 28 February and 20 March 2020. Estimates of the impact of these interventions on SARS-CoV-2 transmission are critical for decision making in this and future outbreaks. We here aim to assess the impact of these NPIs on disease transmission by estimating changes in the basic reproduction number (R0) at national and cantonal levels in relation to the timing of these NPIs. We estimated the time-varying R0 nationally and in eleven cantons by fitting a stochastic transmission model explicitly simulating within-hospital dynamics. We used individual-level data from more than 1000 hospitalised patients in Switzerland and public daily reports of hospitalisations and deaths. We estimated the national R0 to be 2.8 (95% confidence interval 2.1–3.8) at the beginning of the epidemic. Starting from around 7 March, we found a strong reduction in time-varying R0 with a 86% median decrease (95% quantile range [QR] 79–90%) to a value of 0.40 (95% QR 0.3–0.58) in the period of 29 March to 5 April. At the cantonal level, R0 decreased over the course of the epidemic between 53% and 92%. Reductions in time-varying R0 were synchronous with changes in mobility patterns as estimated through smartphone activity, which started before the official implementation of NPIs. We inferred that most of the reduction of transmission is attributable to behavioural changes as opposed to natural immunity, the latter accounting for only about 4% of the total reduction in effective transmission. As Switzerland considers relaxing some of the restrictions of social mixing, current estimates of time-varying R0 well below one are promising. However, as of 24 April 2020, at least 96% (95% QR 95.7–96.4%) of the Swiss population remains susceptible to SARS-CoV-2. These results warrant a cautious relaxation of social distance practices and close monitoring of changes in both the basic and effective reproduction numbers.

Optimal control of the spatial allocation of COVID-19 vaccines: Italy as a case study.

While campaigns of vaccination against SARS-CoV-2 are underway across the world, communities face the challenge of a fair and effective distribution of a limited supply of doses. Current vaccine allocation strategies are based on criteria such as age or risk. In the light of strong spatial heterogeneities in disease history and transmission, we explore spatial allocation strategies as a complement to existing approaches. Given the practical constraints and complex epidemiological dynamics, designing effective vaccination strategies at a country scale is an intricate task. We propose a novel optimal control framework to derive the best possible vaccine allocation for given disease transmission projections and constraints on vaccine supply and distribution logistics. As a proof-of-concept, we couple our framework with an existing spatially explicit compartmental COVID-19 model tailored to the Italian geographic and epidemiological context. We optimize the vaccine allocation on scenarios of unfolding disease transmission across the 107 provinces of Italy, from January to April 2021. For each scenario, the optimal solution significantly outperforms alternative strategies that prioritize provinces based on incidence, population distribution, or prevalence of susceptibles. Our results suggest that the complex interplay between the mobility network and the spatial heterogeneities implies highly non-trivial prioritization strategies for effective vaccination campaigns. Our work demonstrates the potential of optimal control for complex and heterogeneous epidemiological landscapes at country, and possibly global, scales.

Characterizing human mobility patterns in rural settings of sub-Saharan Africa.

Human mobility is a core component of human behavior and its quantification is critical for understanding its impact on infectious disease transmission, traffic forecasting, access to resources and care, intervention strategies, and migratory flows. When mobility data are limited, spatial interaction models have been widely used to estimate human travel, but have not been extensively validated in low- and middle-income settings. Geographic, sociodemographic, and infrastructure differences may impact the ability for models to capture these patterns, particularly in rural settings. Here, we analyzed mobility patterns inferred from mobile phone data in four Sub-Saharan African countries to investigate the ability for variants on gravity and radiation models to estimate travel. Adjusting the gravity model such that parameters were fit to different trip types, including travel between more or less populated areas and/or different regions, improved model fit in all four countries. This suggests that alternative models may be more useful in these settings and better able to capture the range of mobility patterns observed.

The epidemicity index of recurrent SARS-CoV-2 infections.

Several indices can predict the long-term fate of emerging infectious diseases and the effect of their containment measures, including a variety of reproduction numbers (e.g. [Formula: see text]). Other indices evaluate the potential for transient increases of epidemics eventually doomed to disappearance, based on generalized reactivity analysis. They identify conditions for perturbations to a stable disease-free equilibrium ([Formula: see text]) to grow, possibly causing significant damage. Here, we introduce the epidemicity index e0, a threshold-type indicator: if e0 > 0, initial foci may cause infection peaks even if [Formula: see text]. Therefore, effective containment measures should achieve a negative epidemicity index. We use spatially explicit models to rank containment measures for projected evolutions of the ongoing pandemic in Italy. There, we show that, while the effective reproduction number was below one for a sizable timespan, epidemicity remained positive, allowing recurrent infection flare-ups well before the major epidemic rebounding observed in the fall.

Range of reproduction number estimates for COVID-19 spread.

To monitor local and global COVID-19 outbreaks, and to plan containment measures, accessible and comprehensible decision-making tools need to be based on the growth rates of new confirmed infections, hospitalization or case fatality rates. Growth rates of new cases form the empirical basis for estimates of a variety of reproduction numbers, dimensionless numbers whose value, when larger than unity, describes surging infections and generally worsening epidemiological conditions. Typically, these determinations rely on noisy or incomplete data gained over limited periods of time, and on many parameters to estimate. This paper examines how estimates from data and models of time-evolving reproduction numbers of national COVID-19 infection spread change by using different techniques and assumptions. Given the importance acquired by reproduction numbers as diagnostic tools, assessing their range of possible variations obtainable from the same epidemiological data is relevant. We compute control reproduction numbers from Swiss and Italian COVID-19 time series adopting both data convolution (renewal equation) and a SEIR-type model. Within these two paradigms we run a comparative analysis of the possible inferences obtained through approximations of the distributions typically used to describe serial intervals, generation, latency and incubation times, and the delays between onset of symptoms and notification. Our results suggest that estimates of reproduction numbers under these different assumptions may show significant temporal differences, while the actual variability range of computed values is rather small.

The geography of COVID-19 spread in Italy and implications for the relaxation of confinement measures.

The pressing need to restart socioeconomic activities locked-down to control the spread of SARS-CoV-2 in Italy must be coupled with effective methodologies to selectively relax containment measures. Here we employ a spatially explicit model, properly attentive to the role of inapparent infections, capable of: estimating the expected unfolding of the outbreak under continuous lockdown (baseline trajectory); assessing deviations from the baseline, should lockdown relaxations result in increased disease transmission; calculating the isolation effort required to prevent a resurgence of the outbreak. A 40% increase in effective transmission would yield a rebound of infections. A control effort capable of isolating daily  ~5.5% of the exposed and highly infectious individuals proves necessary to maintain the epidemic curve onto the decreasing baseline trajectory. We finally provide an ex-post assessment based on the epidemiological data that became available after the initial analysis and estimate the actual disease transmission that occurred after weakening the lockdown.

Biological fluid dynamics of airborne COVID-19 infection.

ABSTRACT: We review the state of knowledge on the bio-fluid dynamic mechanisms involved in the transmission of the infection from SARS-CoV-2. The relevance of the subject stems from the key role of airborne virus transmission by viral particles released by an infected person via coughing, sneezing, speaking or simply breathing. Speech droplets generated by asymptomatic disease carriers are also considered for their viral load and potential for infection. Proper understanding of the mechanics of the complex processes whereby the two-phase flow emitted by an infected individual disperses into the environment would allow us to infer from first principles the practical rules to be imposed on social distancing and on the use of facial and eye protection, which to date have been adopted on a rather empirical basis. These measures need compelling scientific validation. A deeper understanding of the relevant biological fluid dynamics would also allow us to evaluate the contrasting effects of natural or forced ventilation of environments on the transmission of contagion: the risk decreases as the viral load is diluted by mixing effects but contagion is potentially allowed to reach larger distances from the infected source. To that end, our survey supports the view that a formal assessment of a number of open problems is needed. They are outlined in the discussion.

Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures.

The spread of coronavirus disease 2019 (COVID-19) in Italy prompted drastic measures for transmission containment. We examine the effects of these interventions, based on modeling of the unfolding epidemic. We test modeling options of the spatially explicit type, suggested by the wave of infections spreading from the initial foci to the rest of Italy. We estimate parameters of a metacommunity Susceptible-Exposed-Infected-Recovered (SEIR)-like transmission model that includes a network of 107 provinces connected by mobility at high resolution, and the critical contribution of presymptomatic and asymptomatic transmission. We estimate a generalized reproduction number ([Formula: see text] = 3.60 [3.49 to 3.84]), the spectral radius of a suitable next-generation matrix that measures the potential spread in the absence of containment interventions. The model includes the implementation of progressive restrictions after the first case confirmed in Italy (February 21, 2020) and runs until March 25, 2020. We account for uncertainty in epidemiological reporting, and time dependence of human mobility matrices and awareness-dependent exposure probabilities. We draw scenarios of different containment measures and their impact. Results suggest that the sequence of restrictions posed to mobility and human-to-human interactions have reduced transmission by 45% (42 to 49%). Averted hospitalizations are measured by running scenarios obtained by selectively relaxing the imposed restrictions and total about 200,000 individuals (as of March 25, 2020). Although a number of assumptions need to be reexamined, like age structure in social mixing patterns and in the distribution of mobility, hospitalization, and fatality, we conclude that verifiable evidence exists to support the planning of emergency measures.